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The detail structure of the wave function is analyzed at various refinement levels using the methods 
of wavelet analysis. The eigenvalue problem of a model system is solved in granular Hilbert spaces, 
and the trajectory of the eigenstates is traced in terms of the resolution. An adaptive method is 
developed for identifying the fine structure localization regions, where further refinement of the wave 
function is necessary. 



I. INTRODUCTION 

Multiresolution or wavelet analysis (MRA) 0,0 is in- 
creasingly used in advanced data compression and image 
coding algorithms 0, and several promising approaches 
applied this technique in describing electronic structures. 
Arias el al. elaborated a Kohn-Sham equation based 
density functional (DFT) method, Goedecker and Ivanov 
have also successfully used wavelets for the solution 
of Poisson's equation. There are also attempts to extend 
wavelet calculations to three dimensional structures in a 
computationally manageable manner 

The basic principle of the application of MRA is the 
recognition that the details of the electronic structure are 
not distributed equally over different parts of the system. 
In data compression the details of the picture are system- 
atically included by consecutive refinements of the image 
only in those spatial regions where it is necessary. A 
similar approach to the electronic wave function suggests 
that the fine details of the distribution are concentrated 
around the nuclear cusps and singularities of the two- 
electron density matrix. It has been shown, that the 
surroundings of a molecule can be described at a rather 
rough resolution level ■ We have also demonstrated Q , 
that the extremely hardly representable electron-electron 
cusp singularity of the two-electron density operator can 
be easily reproduced by the method of multiresolution 
analysis. 

The applicability of these approaches strongly depends 
on the necessary depth of refinement. Clearly, a uniform 
refinement of the distributions leads to an exponential 
increase of the computational efforts with growing reso- 
lution level. The region where the refinements are applied 
has to be confined, in order to achieve computationally 
manageable algorithms. In this contribution we study 



some questions concerning several aspects raised by the 
above considerations. We analyze the detail structure 
of the wave function, we explore how the spatial regions 
where further refinement is necessary are to be identi- 
fied, and how can the extension of the eigenvalue problem 
to new subspaces be restricted. We trace how the pro- 
jection of the wave function to different resolution level 
subspaces is changing in the course of the consecutive 
refinement steps. The performance of the adaptive re- 
finement algorithm is tested by numerical comparisons 
to the exact analytical solution of the harmonic oscilla- 
tor. 



II. RESOLUTION STRUCTURE OF THE 
HILBERT SPACE 

Multiresolution analysis theory ensures that an arbi- 
trary clement of the Hilbert space € H can be exactly 
decomposed into orthogonal components as 

oo 

*( x ) = X/ ct s m( x ) + X! X! dme w ™d x )- (!) 

££Z m=0 £eZ 

Here, the basis functions s^{x) = sq(x — £) are equidis- 
tant shifts of the "mother scaling function" so over a 
coarse grid of spacing 1. Refinements to the first ap- 
proximation described by the first summation in Q 
are introduced by the m = level wavelets wqi(x) = 
wq(x — £). The "mother wavelet" wo is orthogonal to 
s , and generally, (s i\sa k } = &tk, (soe\w k) = and 
(wQi\wak) — 5 ik- Further refinements arc expanded by 
the basis vectors w m e(x) — 2 m / 2 wo(2 m x — £) for resolu- 
tion levels m — 1,2,..., with the orthogonality relations 
{sot\w mk ) = and (w ml \w jk ) = 8 m j 5 tk . The subspaces 
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Vb = span{s £|^ G Z}, W m = span{w m ^|^ S Z} consti- 
tute a complete decomposition of the Hilbcrt space 

H = V Q ®W ®W 1 ®--- , (2) 

i.e., expansion is exact. 

We define the resolution structure of the wave function 
by the series of components 

Po* = s ° e 

Qo* = ^2 d oe w oe 

fez 

Qi* = J^dwww, (3) 



where pj and Q m are orthogonal projection operators 
to subspaces Vb and W m , respectively. As according to 
Parseval's equality 

oo 

iip *ii 2 + J2 \\Q™n 2 = i, (4) 

m=0 

we can measure the mth level complexity of the wave 
function by the number ||Q m ^|| , which characterizes 
how important is the detail space W m in expanding 5*. 

The decision, to include or omit W m is crucial in de- 
veloping useful algorithms for the following reasons. The 
basis functions of W m are "sitting" on an equidistant grid 
with a grid length of ~ 2~ m . If a function constrained to 
a finite domain of the space is expanded, the number of 
basis functions in W m is increasing like ~ 2 mD , where D 
is the dimension of the system. This exponential "explo- 
sion" makes the direct application of ||TJ unacceptable. 
Detail spaces with negligible ||<2 m \P|| 2 (lower than a pre- 
defined threshold) can be completely ignored. The expe- 
rience 0, @ of finding details in a constrained region of 
the space leads, however, to the conclusion, that even if 

\\Q m n 2 =J2 d ™t ( 5 ) 

is significant, only very few terms in the summation 
contribute essentially to its value. This recognition 
would help to avoid the exponential explosion men- 
tioned above, by using the restricted detail space W m = 
span{w TO ^| where d m i is significant}. 

The outlined strategy is, however, useless, if the choice 
of the significant d m e would be based on a prior calcula- 
tion of all coefficients, and testing how many d m e are nec- 
essary to fulfill (JSJ with a good approximation. Clearly, 
a predictive method is needed. 

A. Decomposition of exact wave functions 

We will illustrate the above consideration with a sim- 
ple exactly solvable example in D = 1 dimension. The 



case of D = 3 is expected to have similar behavior with 
correspondingly greater numbers of basis functions. The 
ground state and some excited states of the standard ex- 
ample of the ID harmonic oscillator with the Hamilto- 
nian 

are analyzed. One reason for choosing this system is 
that the matrix elements (soi\H\s k) , (s i\H\wjk) and 
(w m i\H\wjk) can be calculated exactly for |J?J|, making 
possible to avoid inaccuracies not connected to Hilbcrt 
space constraints. The grid length of the scaling function 
subspace Vq was set to 1 a.u. The expansion coefficients 
were calculated as 

ct = (s(m|*) and d m i = (7) 

where *f> stands for the exact ground or excited state 
wave function. The Daubechies-6 [jj set was chosen 
for the multiresolution basis set {sqi, w m e}. The ele- 
ments of this set have finite support (the mother scal- 
ing function is zero outside of the interval [0,5)), and 
both the scaling functions and the wavelets are differ- 
entiable. The scalar products in J7J) were numerically 
calculated on a grid of 2~ 15 a.u. density. Fig. shows 
the values ||Q m Wj|| 2 for various excitations i — 0, . . . , 5 
with to = 1, whereas the projections to the scaling func- 
tion subspace are given as follows: IjPo^oll 2 = 0.9972, 
HP,*!!! 2 = 0.9822, ||P * 2 || 2 = 0.9561, ||Po* 3 || 2 - 
0.8637, ||Po* 4 || 2 = 0.9362 and ||P * 5 || 2 = 0.5582. It 

IIQm**|| 2 
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FIG. 1: The amplitudes of the projections of the exact eigen- 
functions of the Hamiltonian (JSJ with u = 1 to the detail 
spaces W m for m = 0, . . . , 6. The sign stands for the ground 
state i — 0, whereas the signs x, *, o, + and □ denote the ex- 
citations i = l,2, 3, 4, 5, respectively. Atomic units are used. 

is clear, that the bulk of the states is described at level 
to = 0, and the contributions of the detail spaces disap- 
pear with increasing to. It is important to realize, that 
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even if the number of basis functions in subspace W m in- 
creases as ~ 2 m , the total contribution of the subspace to 
the wave function exponentially decreases. This results 
coincides with our previous statement that many- 
electron density operators can not contain arbitrary fine 
(nor rough) details. 



B. Fine structure localization 

Besides a "vertical" truncation of the Hilbert space 
over a sufficiently large resolution level M, there is a 
possibility of reducing the size of the subspaces W m in 
a "horizontal" truncation process, by decimating the ba- 
sis functions w m g which belong to those spatial regions 
where the wave function does not contain fine details. 

In order to study the extent of this fine structure 
localization, we have also examined, how many coeffi- 
cients a and djnt are essential in the norms ||Po^i|| 2 an d 
||<2m^i|| 2 ! respectively. The projections of the wave func- 
tion to the restricted subspaces Vq and W m are defined 
by the projectors Pq and Q m . 

After numerically calculating the scalar products JJJ , a 
threshold value 77, close to 1 was chosen. We counted the 
number #Vb and jfW m of most significant coefficients for 
which the inequalities 



l|p *iir- \\Po*i\f < i-v 

\Qm%\\ 2 ~ \\Qm%\\ 2 < 1-r? 



and 



(8) 



hold. This criterion allows each subspace W m to intro- 
duce an error 1 — 77 uniformly, and for those m where 
II Qm^i || 2 < 1 — 77 the complete subspace can be omit- 
ted, thus W m = 0. Fig. |2 summarizes the results for the 
ground and a selected excited state as a function of rj 
for various resolution levels. On the horizontal axes the 
number of digit 9 in r\ is shown, i.e., the value 1 corre- 
sponds to rj = 0.9, whereas 5 corresponds to 77 = 0.99999, 
etc. It is clearly seen, that ground and excited states be- 
have in a very similar manner. As the logarithmic plots 
show, in the rj — > 1 limit the number of basis functions 
in the subspaces W m increases as ~ 2 m . On the other 
hand, if a finite, but even relatively high precision, like 
rj = 0.9999999999 is demanded, many of the high res- 
olution detail spaces drop out, and even at lower reso- 
lution m, only an acceptable amount of basis functions 
contribute to the truncated subspaces. 



III. TRAJECTORY AND LIMITING BEHAVIOR 
OF THE EIGENFUNCTIONS IN 
COARSE-GRAINED HILBERT SPACES 

In applied calculations the granularity M of the Hilbert 
space (i.e., the maximum resolution m in the summations 
of and can not be arbitrary large. Instead, for a 
finite granularity M the eigenvalue problem is solved in 



the restricted Hilbert space 

H [M] = V ®W ®W 1 ®---®W M -i, (9) 

which assumption seems to be warranted in the light of 
the above results. For M = 0, we define H [0] = Vq. The 
solution of the eigenvalue problem of H in the coarse- 
grained Hilbert space H\ M ^ leads to the ground and ex- 
cited states <J/t M ' (i = 0, 1, 2, . . .). In this case we lose fine 
details, thus the natural question arises, how the trajec- 
tory ^f 1 — > *f — * \ff[ 2 1 — ► • ■ • approaches the limit ^j. 

We have determined the eigenfunctions of the Hamil- 
tonian JBJ with u = 1 for M = 0, . . . , 5. In order to keep 
the number of basis functions finite, the spatial exten- 
sion of the system has to be limited. We have chosen 
\x\ < x max = 10 in ©. The error made by this assump- 
tion can be estimated by the omitted part of the norm 
square of the ground state 



1 - 



|*o| 2 dx = l- erf(Vw£ 



which is less than the accuracy of the double precision 
numerical representation. 

The four-index matrix elements of H can be reduced to 
two-index quantities by a simple variable transformation. 
The kinetic energy matrix elements are calculated using 
the fact that the momentum operator is self-adjoint, thus 
for mi < TO2 

(soiA - V 2 /2 \s £ 2 ) = ^ 1 {s' m \s' til2 _ ll ), 

{ SUl \~W 2 /2\w mA ) = 2- 1 ( S , Q0 \w' m2L ), (10) 

(«Wi|-V a /2|«W a ) = 1 2mi ~ l « Q W m2 - miL ), 

with L = £ 2 - 2 m2 ~ mi £i. The potential energy matrix 
elements are 

(s Qei \uj 2 x 2 /2\s i 2 ) = aj 2 2^ 1 ((s oa \x 2 \s i 2 -i 1 ) 

+2£i(s o\x\soi 2 -i 1 ) 
+tl{s o\s z 2 -e 1 )), 

{s Ml \u 2 X 2 /2 \w m2 £ 2 ) = OJ 2 2~ 1 ((s o\x 2 \w m2L ) 

+2£i(s 00 \x\w m2 1) 

+ei(s o\w m2 £>) » (11) 

{w mitl \u 2 x 2 /2\w m2t2 ) = u 2 2- 2m '- x - 

— mi LI 
+2£i ( Wqq \ x\Wm 2 — mi l) 

+£ 2 (woo|w m2 - mi L)). 

At any resolution level, the x 2 and x operators can 
be exactly expanded as linear combinations of scal- 
ing functions Wavelets are also finite linear com- 
binations of scaling functions of the succeeding res- 
olution level. Consequently, the remaining part of 
the calculation contains the evaluation of integrals 
of type J s f _ ll x 1 (x)s^ 2 x 2 (x)s fl3 x 3 (x)dx with s^x(x) = 
2 >1 / 2 So(2 il x — A). There exist special algorithms for this 
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FIG. 2: The number of the significant basis functions in various subspaces required to reproduce the 1 1 i-fe ^ i 1 1 2 and ||Q m ^/i|| 2 to 
a given precision. The horizontal axis shows the number of digit 9 in the threshold value rj. The sign + stands for the restricted 
scaling function subspace Vo, whereas the signs o, □ , 0, V, A, ~k and * denote the detail spaces W m , with m = 0, 1, 2, 3, 4, 5, 6, 
respectively. The results for the ground state (i = 0) are plotted in the left column both in linear and log scale, while the right 
column corresponds to the third excited state. 
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task, and for exactly determining the integrals of prod- 
ucts of scaling function derivatives [13 . Finally, the 
Hamiltonian matrix was diagonalized by standard sub- 
routines. 

At each refinement level (granularity) M the diago- 
nalization of the Hamiltonian matrix ifL^M] restricted 
to the appropriate subspace Ti\ M ^ leads to the eigenvec- 
tor components c^g and S^i ■> which are the expansion 
coefficients of 

M-i 

= E4** 1 •«(*) + E «ww (12) 

for the approximate eigenstates. As the basis set of 
n [M] = n [M-i] w M -i is an extension to that of 
j^lM—i] ^ expansion coefficients resulted from indepen- 
dent diagonalizations of H\ H [m-i\ andif|^[M] can change 
in the refinement step M—l — > M. Describing the trajec- 
tory of "I'M in 7i in terms of the granularity, would need 
to trace all the coefficients cjjj; and d^, this is how- 
ever, hard to demonstrate. Instead, we have collected 
the ground state projections to a given detail space, and 
presented them in Tab. [I] The excellent convergence 
can be easily realized. After few refinement steps the 
coefficients stabilize. 

For higher excitations the overall trend is similar to 
that illustrated in Tab. U for the ground state, with the 
remark, that the contributions of finer detail spaces are 
more and more emphasized for increasing excitation lev- 
els. As an illustration, we recall the projections of the 
i = 3 excited state to the wavelet subspaces at resolution 
level M = 5: ||P *3 51 || 2 = 0.86367, ||Qo* 3 51 || 2 = 0.13284, 
IIQi*! 51 !! 2 = 3.4188 x 1(T 3 , ||Q2^ 5l || 2 = 6.5664 x 1(T 5 , 
||Q 3 il4 5l |j2 = 1.0810 x 1CT 6 and \\Q^f \\ 2 = 1.7114 x 

io- 8 . 

For completeness, we have also given the errors of the 
ground and excited state energy approximations in the 
Hilbert spaces ft[ M l in Fig. Exponential convergence 
in terms of the granularity level M can be clearly identi- 
fied. 



IV. AN ADAPTIVE METHOD FOR 
DETERMINING FINE STRUCTURE 
LOCALIZATION REGIONS 

In realistic cases, of course, exact solutions of the eigen- 
value problem are not known. We can suppose, how- 
ever, that quantum mechanical systems behave similarly 
to the simple exact case studied above. We expect, that 
high resolution wavelets can completely be omitted, and 
even at lower resolutions only a fractional part of them 
contribute essentially to the expansion. In the following 
considerations we will apply a consistent notation for the 
two distinct truncation schemes. Upper index [M] de- 
notes the "vertical" cutoff of the Hilbert space above the 
granularity M, whereas tilde is applied to indicate the 
"horizontal" truncation in a given detail space W m . 




10 ' 1 1 1 1 1 

1 2 3 4 5 

M 

FIG. 3: The errors of the approximate energies E\ M ^ deter- 
mined in the Hilbert space 7i' M '. The sign stands for the 
ground state i — 0, whereas the signs x, *, o, + and □ denote 
the excitations i = 1,2, 3, 4, 5, respectively. Atomic units are 
used. 



Our goal is to develop a method for determining the 
significant coefficients in JIJ, without calculating all the 
negligible ones. This requires predicting the set of im- 
portant fine level coefficients for M + 1, supposing that 
the approximate expansion 

M-i 

= e 4 M1 + E E ^ 

iev m=o eew m 

(13) 

of the wave function is already known. Here, we have 
used the loose notation I S W m , indicating the fact, that 
the summation index I is constrained to the wavelets w m i 
which belong to the restricted subspace W m . In the first 
summation Vq denotes the restricted subspace of those 
scaling functions which are localized in the finite domain 
of the space, which is occupied by the system. 

The iterative extension vfrI A/ ] — > \j/[ M +!] ca n be sep- 
arated to two independent problems. The first step is 
determining the significant part Wm of the Mth level de- 
tail space, knowing the approximation \E»[ M 1 . In the other 
step, the numerical determination of the expansion coef- 
ficients (1m i in Wm should be carried out. In this work, 
we will concentrate on the solution of the first problem. 
In order to distinguish the errors emerging from different 
sources, the second part of the calculation is carried out 
by an exact diagonalization of the Hamiltonian jj[ M+1 ] 
restricted to the subspace H [M+1] =V ®W ®---® W M - 
In optimized numerical works sophisticated methods for 
solving the eigenproblem of sparse matrices can be ap- 
plied. 

As \E»[ M ] is the eigenfunction of the restricted Hamil- 



TABLE I: The trajectory of the ground state wave function in the detail spaces Vo, Wo, ■ ■ ■ , the contribution of the detail spaces to the Mth coarse-grained 
approximations of the ground state eigenfunction. The last column shows the sum of the first five terms in Parseval's equation |J3J. In finite subspaces 7i' ', . . . , Tv 6 ' 
the relation fulfills exactly, whereas in the infinitely fine Hilbert space the error is of order 10 -12 . 
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tonian if l M I = H\^ [M] , it satisfies the equation 

(fiM_£[«l)$[«l =0 . (14) 

This condition does not hold, if cither or 
are replaced by their exact counterparts, and the qual- 
ity of the Mth resolution level approximate solution of 
the eigenvalue problem can be characterized by the error 
function 

(H - EM)y[M] ^ o. (15 ) 

Based on the methods developed in 0, , we suggest 
the following criterion for selecting the important basis 
functions of Wm- The error function (1 1511 will be mea- 
sured in the extended Hilbert space 7YJ M+1 ] (instead of 
the full Hilbert space TC). Considering, that according 
to 11411 the expansion coefficients of the error function 
(w me \(H - EM^Wl) = if w ml e U [M \ the magni- 
tude of the deviation from the exact eigenvalue problem 
can be characterized by the magnitude of the expansion 
coefficients in the detail space Wm 

r Mi = \{w M( \H-E^\^)\ 2 . (16) 

The wavelet wmi is considered to be "important" if Tmi 
is larger then a given threshold. This test, however, can 
not be carried out in a reasonable time for all the basis 
functions in Wm, as this would require an exponential 
amount of work. Here, we apply an adaptive approach. 
We test the wavelets wmi only for those values of £, where 
the support of wmi overlaps with the previously found 
"important region" of Wm-i- 

To be more specific, we apply the following procedure. 

1. We select a preliminary set Wm of wavelets id mi 
for which the condition 

support iwMi) H support (Wm-i) $ 

holds. 

2. For the selected wavelets the values of (|16|) are cal- 
culated and the wavelets are sorted according to 
descending order of Tmi- 

3. The first most important wavelets constituting the 
restricted subspace Wm are kept until the condition 

^ r M i - 51 r M l < 1 - V, ( 17 ) 

eew M iew M 

similar to (JSJ is fulfilled. Depending on the required 
precision r\ the set Wm is usually much smaller then 
the candidate set Wm- 

The adaptive fine structure localization method reduces 
the size of the detail spaces considerably, for example, 
in the case of nine-9 precision (i] = 0.999999999) and 



the ground state wave function, dim Wq = 646, whereas 
dim We = 296. 

The results obtained for the wave function with various 
precision requirements ij are summarized in Tables Inland 
IIIII The wave function vpI M l calculated in the "horizon- 
tally" truncated Hilbert space is expected to significantly 
deviate from the eigenstate ^[ M l received without trun- 
cation if the required precision is low. On the other hand, 
in the 77 — >• 1 limit \jH M l approximates $[ M 1 very well. 

Tab. HJ contains the deviations A||P *o MI || 2 = 
||P *] M1 || 2 ll^o^H 2 and A||Q m $Hl 2 = 
\\Qm^ l M] \\ 2 - \\Qm^ l o M] \\ 2 for a rough precision 
r\ = 0.99. It can be seen, that the error of the wave 
function saturates with increasing resolution level M. 
The value of the error is significantly larger than the 
error caused by the "vertical" truncation of the Hilbert 
space. 

Tab. lIIll shows similar data for a moderate precision re- 
quirement 77 = 0.999999. In this case the wave function 
<f[ M ], obtained using the adaptive fine structure localiza- 
tion method, gives an excellent approximation to \I/[ M 1. 
The error introduced by the "horizontal" truncation is 
less, than that of the "vertical" one, up to the resolution 
level M = 5. 

In case of excited states, the eigenstates of the eigen- 
value problem in granular Hilbert spaces with a given 
granularity M are less accurate than the ground state. 
The deviations due to "horizontal" truncation of the de- 
tail spaces, however, only slightly exceed the ground state 
values. Consequently, for a given M, the precision r\ 
which was appropriate for ground state calculations will 
certainly be applicable in excited state calculations, as 
well. 

The errors of the approximate energies calculated for 
the ground and excited state i — 3 are plotted in Fig. 0] 
for various precision requirements 7/. For lower precision 
the magnitude of the error saturates at various resolu- 
tion levels M, depending on the value of t). This fact, 
together with the observation concerning the error of the 
wave function, emphasizes the importance of choosing 
matching values for the precision requirement 77 and the 
granularity level M. 



V. CONCLUSIONS 

We have shown, that the fine structure of the Hilbert 
space vanishes exponentially in realistic wave functions. 
Realizing this fact, one can obtain excellent approxima- 
tions of the eigenfunctions by omitting the fine resolution 
detail spaces and solving the eigenvalue problem in gran- 
ular Hilbert spaces. We have traced the trajectories of 
the approximate wave functions of an exactly solvable 
model system, and concluded, that they approach the 
exact wave function exponentially fast. A similar state- 
ment is valid for the approximate energy values. 

A further possibility for reducing the computational 
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TABLE II: The deviation of the trajectory of ^ M ' from that of \I> M ' due to fine structure localization with r\ 



0.99 in the 



detail spaces Vq, Wo, 



,W 4 



M 



[M] II 2 



AiiQ §nr 



[M] 1 1 2 



[M] 1 1 2 



0.41770 x 10~ 5 
3.22557 x 10 -5 
3.14686 x 10 -6 
2.76469 x 10~ 5 
2.75634 x 10~ 5 



0.41770 x 10" 
3.22325 x 10" 
3.17587 x 10" 
2.82869 x 10" 
2.81104 x 10" 



-0.23237 x 10~ 7 

2.50837 x 10~ 7 3.92755 x 10~ 8 

5.95967 x 10~ 7 4.31415 x 10~ 8 

5.11324 x 10~ 7 3.44685 x 10~ 8 



0.89745 x 10" 
1.19115 x 10" 



2.901 x 10" 



TABLE III: The deviation of the trajectory of from that of * due to fine structure localization with r) = 0.999999 in 
the detail spaces Vb, Wo, • • • , W4. 
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AiiPo*nr 


AiiQ *nr 


AiiQi*nr 


A||Q 2 * M] | 
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A|iQ 3 *nr A|iQ 4 *nr 
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-0.00209 x 10 -9 


0.00209 x 10 -9 
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0.21649 x 10~ 9 


-0.21870 x 10~ 9 


0.0220 x 10~ 10 
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-1.02284 x 10 -9 


1.04319 x 10~ 9 


-0.1999 x 10 -10 


-0.036 x 10~ 


-11 




4 


4.15413 x 10~ 9 


-4.42360 x 10~ 9 


2.9803 x 10~ 10 


-2.845 x 10" 


-11 


-1.0 x 10~ 13 


5 


3.26705 x 10 -9 


-3.50721 x 10~ 9 


2.6070 x 10~ 10 


-2.072 x 10" 


-11 


1.8 x 10~ 13 -10~ 14 
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FIG. 4: The errors of the approximate energies s] M ' determined in the Hilbert space 7i' M ' for the ground state i = and for 
an excited state i = 3. The signs x, □, 0, *, A and o stand for the threshold values 77 = 0.9, 0.99, 0.999, 0.9999, 0.99999 and 
0.999999, respectively. Atomic units are used. 



complexity of the calculations, that higher resolution 
wavelets are included in the basis set only in those spa- 
tial regions where the fine structure of the wave function 
requires it. For our model system we found, that the 
number of significant basis functions is considerably less 
than the dimension of the detail subspaces included in 
full calculations. The a priori selection of the significant 
basis functions needs a predictive algorithm. 



Based on the above concept, we have developed an 
adaptive method for selecting the essential basis func- 
tions using the fine structure localization technique. Cal- 
culations in truncated Hilbert spaces, restricted in such 
manner lead to sufficiently precise wave functions and 
eigenenergies, even in the case of moderately strict basis 
function selection criterion. 
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